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ABSTRACT 


Transverse MHD bending waves are considered in an isothermal 
and compressible two-dimensional current sheet of finite 
thickness in which the magnetic field changes direction and 
strength. The general form of the wave equation is obtained. It 
is shown that rotation of the magnetic field across the current 
sheet prevents the existence of singular points so that 
continuous spectrum solutions and the concomitant wave decay 
disappear. Instead, normal modes exist and closed integral 
solutions for arbitrary current sheet structure are found. 

The results are discussed in terms of small-scale waves on 
the heliospheric current sheet. 


Subject headings: sun:solar wind - hydromagnetics - surface 


waves 
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1 . INTRODUCTION 

Bending waves are an essentially transverse displacement of 
an equilibrium interface or layer in which, to first 
approximation, gradients along the layer are negligible in 
comparison to gradients across the layer. Recently, Bertin and 
Coppi (1985) applied the theory of MHD bending waves to a current 
sheet in which the magnetic field reverses direction. They 
suggested that these waves might provide an explanation for the 
well-known sector structure in the interplanetary magnetic field 
polarity (Wilcox and Ness, 1965) which has often been interpreted 
as a warping of the heliospheric current sheet (HCS) dividing 
opposite magnetic polarity regions in the interplanetary medium 
(Hundhausen, 1972, op . ci t . ; Suess and Hildner, 1985). 

It is now generally accepted that warping of the HCS does 
exist but that it is a direct consequence of the large-scale 
structure of the solar magnetic field as that field is carried 
outward into the interplanetary medium by the solar wind 
(Hoeksema, et al . , 1982, 1983; Hoeksema, 1984). The large-scale 
warping is therefore understandable in the kinematic limit where 
the internal dynamics of the current disk is altogether 
ignorable. A surface plot of the current sheet made using this 
approximation is shown in Figure 1. The plot was made under the 
assumption that the interplanetary field is the projection of a 
"source surface" potential magnetic field model of the corona 
ending at 2.5 solar radii with a constant solar wind velocity of 
400 km/s. The epoch chosen is the time the Giotto spacecraft 
passed by Halley's comet - approximately the same instant that 
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both objects passed through the current sheet. In situ 
verification has shown that these projections are valid to within 
one day at 1 AU. 

Nevertheless, important reasons exist for studying bending 
waves on the HCS. There are small-scale "ripples" on the HCS 
that have so far eluded complete understanding. Consequently, 
the fine-scale structure of the HCS is a topic of active 
investigation (Behannon, et al . , 1981). It is worth noting that 
the analysis by Bertin and Coppi (1985) explicitly invoked a WKB 
approximation that, a priori , made the results far more relevant 
to these small-scale fluctuations in the first place and 
seemingly invalidated their discussion in terms of the large- 
scale sector structure 

Here we re-analyze bending waves in the ideal MHD limit in a 
two-dimensional current sheet of finite thickness - implicitly 
invoking the WKB approximation by ignoring gradients along the 
current sheet. The unperturbed state will be in equilibrium with 
its surroundings, isothermal, and compressible. The waves will 
also be assumed to be isothermal but incompressible; thus, we 
consider the Alfven type of purely transverse waves. These are 
the same assumptions made by Bertin and Coppi (1985). With no 
additional complexity, we are able to consider any two- 
dimensional, equilibrium current sheet in which the magnetic 
field changes direction and/or strength, with a current sheet 
across which the field reverses direction (a "neutral sheet") in 
an arbitrary manner being a special case. 
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Our approach is a generalization of the often-studied 
problem of hydromagnetic surface waves associated with a layer 
between otherwise homogeneous domains in static eguilibrium that 
has been analyzed in many contexts (see Lee and Roberts, 1986, op 
cit . ; Hollweg, 1986; Roberts, 1983). These earlier studies have 
shown that, for a layer in which the magnetic field vector does 
not change direction and the field strength has no local 
extremum, a singular point exists where the wave frequency equals 
the local Alfven frequency. The physical consequence of this is 
the decay of a surface disturbance propagating along the layer. 
The decay results from "mode conversion" (sometimes called "phase 
mixing") of the collective surface disturbance into local 
oscillations within the interface. That is, there is no normal 
mode solution for these wave numbers - a continuous spectrum of 
modes exists except in the limit that the thickness of the layer 
goes to zero. This behavior has been analyzed most accurately 
and clearly by Lee and Roberts (1986). 

Under particular circumstances it is also possible to have 
normal mode solutions for a limited region in (a>,k) space (where 
a) is the circular frequency and k is the wavenumber) in a 
compressible medium (Hopcraft and Smith, 1986). What we show 
here is that in a current sheet in which the field vector 
undergoes rotation in the plane of the sheet, in addition to an 
amplitude change, the region of (w,k) space permitting normal 
mode solutions expands and that continuous mode solutions can 
sometimes disappear altogether, along with the concomitant wave 
decay. This is particularly interesting in reference to the HCS, 
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wherein it is observed that the magnetic field normally reverses 
direction without the amplitude going to zero (Behannon, et al . , 
1981) . 

Our formalism assumes waves with no compressive component in 
an otherwise compressible medium of finite sound speed, as 
opposed to assuming the limit of infinite sound speed as outlined 
by Priest (1982, egn. (4.60)) or assuming an incompressible 
medium. This allows us to reduce the order of the wave equation 
for normal modes and eventually to find, for the first time, a 
closed integral solution for arbitrary current sheet structure 
for those conditions when normal modes exist. Our wave equation 
reduces to a known equation in various limits, giving us 
confidence in our results. 

In the following, we explicitly formulate the bending wave 
problem in section 2, present the MHD bending wave equations in 
section 3, and discuss our results in an astrophysical context 
and in terms of small-scale waves on the HCS in section 4. 


2. FORMULATION OF THE PROBLEM 

In this section, we describe the model and geometry of the 
current sheet, discuss the type of bending waves which can be 
excited in the current sheet, and present the governing 
magnetohydrodynamic (MHD) equations needed to describe these 
waves. We consider the ideal MHD equations obtained by 
neglecting diffusion of the magnetic field, displacement currents 
and electrostatic forces and by making the assumption that the 
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gas pressure is a scalar. We begin with the current sheet model, 
then present the perturbations in the current sheet and the MHD 
equations . 

(a) The Model and Geometry of a Current Sheet 

A current sheet is characterized by changing strength and 
direction of the magnetic field and by the pressure balance 
across the sheet. From the theoretical point of view, both the 
magnetic field and the gas pressure can change across the current 
sheet in an arbitrary way (Bertin and Coppi, 1985), however, for 
a current sheet in static equilibrium the total (gas + magnetic) 
pressure has to be constant. On the other hand, in situ 
observations of the HCS (Behannon et al . 1981), across which the 
field changes polarity, give restrictions on the behaviour of the 
magnetic field across the current sheet? the data show that the 
magnetic field predominantly rotates across the sheet but also 
that some variations of the strength of the field can be seen. 

In the latter case, gas pressure variations across the sheet are 
needed in order to maintain the pressure balance. 

In the approach presented in this paper, we consider an 
isothermal current sheet in static equilibrium. We allow for 
rotation of the magnetic field across the current sheet and for 
arbitrary changes in the field strength; the latter requires 
changes of density across the sheet to account for the pressure 
variations. We impose perturbations that bend the current sheet 
and look for the effects on the behaviour of MHD bending waves 
when gradients of the physical parameters across the sheet are 
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taken into account. Variations of the physical parameters along 
the current sheet are neglected in this approach as they are 
assumed to be on characteristic scales much larger than those 
considered across the current sheet. We introduce a coordinate 
system with the x and y-axis in the plane of the sheet (the y- 
axis is fixed in the direction of the magnetic field on the 
"upper" edge of the sheet) and with the z-axis normal to the 
sheet (see Figure 2a); all the physical parameters that describe 
the medium inside the sheet vary across the sheet and show 
dependence on z alone. 

Having defined the coordinate system, we can, with complete 
generality, describe the rotation and gradient of the magnetic 
field inside the current sheet in the following way: 

S q (z) = B ox (z)x + B Qy (z)y = B q (z)[x sin a(z) + y cos a(z)], 

( 2 . 1 ) 


where a is the angle between the direction of the magnetic field 
and the y-axis (see Figure 2b), and B Q (z) is an arbitrary 
function of z. The static equilibrium is defined by 


7p o - bt ' 7 x 6 o> x = °> 


( 2 . 2 ) 


which can easily be reduced to the pressure balance of the 
current sheet (along the z-axis) with its surroundings 


B 


= P (1> 
*e 


8 IT 


= P (2) 


( B e 2) V 

8 TT 


= const. 


(2.3) 
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where p Q is the gas pressure inside the current sheet, and 

p^l ^ , p^), and ^ are the gas pressure and magnetic field 

outside, on opposite sides of the current sheet. Both latter 

quantities are assumed to be constant in the vicinity of the 

current sheet as we look at wavelength larger than the thickness 

of the current sheet. and B^ ^ need not be in the same 

e e 

direction in the x-y plane. 

Cb) Perturbations of the Current Sheet 

To develop a formalism for describing small-amplitude MHD 
bending waves of arbitrary frequency, we consider perturbations 
of the velocity and the magnetic field normal to the current 
sheet; these perturbations bend the current sheet and are 
transverse MHD bending waves propagating in any directions along 
the current sheet. The amplitude of the density compression is 
negligible for small-amplitude transverse waves and to the first- 
order approximation the driven compressional wave may also be 
neglected. In the approach presented here, we neglect density 
and pressure perturbations. We allow for arbitrary propagation 

. A A 

along the current sheet (the wave vector k = k x x + k^y, see 
Figure 3) and consider the most general form of velocity 
t)(x,y,z,t) and magnetic field §(x,y,z,t) perturbations: 

3(x,y,z,t) = 0 p (z)exp[-i(<ut - k x x - k y y) ] , (2.4) 


and 
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£(x,y,z,t) = $ (z)exp(-i(ojt - k x - k y)], (2.5) 

p x y 

where w is the wave frequency. 

As shown by equations (2.4) and (2.5), the velocity O^tz) 
and the magnetic field ^p(z) perturbations have three components, 
however, only the z-component of the perturbations excite MHD 
bending waves which are purely transverse waves propagating along 
the current sheet. Purely transverse waves that show 
oscillations of the velocity and magnetic field only in the xy- 
plane can also be excited by the x- and y-components of equations 
(2.4) and (2.5); in this case, however, the wave is not a bending 
wave. Note also that this transverse wave cannot be excited by 
the perturbations along the magnetic field as there is no the 
restoring force acting in that direction. 


(c) The MHD Equations 

Our model of the current sheet and the type of the 
perturbations chosen reduce the ideal MHD equations to the 
following linearized form: 


V • 



( 2 . 6 ) 


ft - 47^ " 7 X 6) X ®o + <7 * « 0 > X S] = 0, 

f! - V X <8 X s 0 > = o. 


(2.7) 


( 2 . 8 ) 
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where is the density inside the current sheet. Note that the 
o 

equation of motion is simplified by the static equilibrium 
condition (equation 2.2) and that the perturbations of 6 and 6 
can be calculated from equations (2.7) and (2.8) alone; the 
continuity equation (2.6) and the solenoidal condition (v«§ = 0) 
give additional restrictions on these perturbations. 

The set of MHD equations (2.6) through (2.8) fully describes 
the propagation of purely transverse waves in a compressible 
medium and allows us to generalize the approach developed for 
incompressible ideal MHD (Bertin and Coppi, 1985; Lee and 
Roberts , 198 6 ) . 

In order to consider the effect of gradients inside the 

current sheet on the behaviour of MHD bending waves the wave 

equation for either the velocity U pz (z) or the magnetic field 

B_,(z) perturbation has to be derived and solved. We will 
P *• 

address these problems in the next section. 


3. THE MHD BENDING WAVE EQUATIONS AND THEIR SOLUTIONS 

In the following, we derive the wave equations for MHD 
bending waves, show general solutions and discuss limiting cases. 

(a) The Wave Equations 

In the standard approach, the velocity and magnetic field 
perturbations defined by equations (2.4) and (2.5) are 
substituted into the equation of motion (2.7) and the induction 
equation (2.8). Then, with the help of the continuity equation 



11 


(V»0 = 0), the wave equation for the velocity perturbation U__ 

P P Z 

is obtained in the form 


_d 

dz 


(D 


dU. 


_pz. 


A dz 


k 2 C> U 
A pz 


(3.1) 


where D. = u> 2 - (£•$.) 2 and where = [B x + B y]/(4wp ) ^ 2 . 

A a A ox oy o 

This equation describes transverse MHD bending waves in a fully 
incompressible medium and has been extensively studied by many 
authors (e.g., Tataronis and Grossmann, 1973; Roberts, 1984? 
Bertin and Coppi, 1985; Lee and Roberts, 1986). 

In the approach presented here, we consider a compressible 
medium and combine the continuity equation (2.6) with the 
induction equation. It gives 


dB 


io)B + iU ( £ • § ) - 0 n + U B X = 0, 

px px o pz v dz pz ox r 


dB 


iuB + iU tk.S ) - U f— 3^) + U B A = 0, 
py py o pz'' dz > pz oy r 


(3.2) 

(3.3) 


io>B + iU (£.£ ) = 0, 
pz pz o 


(3.4) 


where A r = (dp o /dz)/p Q . 

To calculate the wave equation for U we eliminate B„„ , 

pz px 

Bpy, and Bp Z in terms of Up Z from equations (.3.2), (3.3), and 
(3.4) using the x and y components of the equation of motion. 
After some algebra, we get 
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iwB 


px,py 


■[t 


dB k ’V? ’ 

If* 2 *) . _2EiJ^*_,t.« n)Xl 


to 


+ ( 


U 2 B - (£.3 )k V* 

P - o X 4 y ----) X B 2 1 U , 

.2 „2 ; r oj pz 


u 2 B 2 - (£.3 ) 


pz' (3.5) 


where X b = (dB 0 /dz)/B Q . 


Having obtained B nv and B_„, we may eliminate them from the 

y A py 

z-component of the equation of motion and obtain the wave 

equation for U „ in the form 
pz 


°A U 


pz 





0, 

(3.6) 


2 2 

where again D A = u - ( £ A ) ^ is the dispersion relation for 
Alfven waves. To further reduce this equation, we use the 
pressure balance across the current sheet (equation 2.3) and find 
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where V c (= p /p ) is the sound speed. Thus, the MHD bending 

o O O 

wave equation for U_._ can finally be written as 


o) 2 D. U 
A pz 


d_ 

dz 


[ ’o V a (d a - 


“M’Vpz 1 = °' 


(3.8) 


where b a = V A / V S - 

Following the same procedure, one may obtain the MHD bending 
wave equation for Bp Z in the form 
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<i) 


Vpz + [" 5 - S c>»o V A (D A - “ 2 *l>»b B pzl ■ 


(3.9) 


which is slightly different than for U pz . It is important to 
know the equation for B pz and its solution as the magnetic field 
can usually be determined more easily than the velocity from 
observational data. 

Both wave equations presented above may be simplified to the 
well-known dispersion relation for the Alfven waves when the 
gradient of magnetic field disappears (x^ = 0) and no layer 
exists. The general solutions can also be easily found in 
analytical forms (see next subsection), however, in order to 
calculate variations of the velocity or magnetic field 
perturbations across the current sheet the explicit form for the 
magnetic field has to be given (see subsection 3e for detailed 
discussion). The rotation of the magnetic field (being forceless 
in the equation of motion) does not play as important role in 
equations (3.8) and (3.9) as gradients of density and magnetic 
field; the latter effects introduce forces into the equation of 
motion and determine the wave behavior. 

Comparison of the wave equations (3.8) and (3.9) to equation 
(3.1) shows that now we are dealing with first order differential 
equations instead of second. Note also that the wave equations 
(3.8) and (3.9) are valid for a compressible medium and have a 
form that is simple to deal with. 
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(b) Singular Points 

In the wave equations for MHD bending waves (equations 3.8 
and 3.9) the derivative with respect to z can be explicitly 
calculated and the wave equations can be written in the following 
forms : 


dU 




pz _ 


r 9 d 6 dD. 

= - [8,p o u 2 D A + D^-gf + & 0 -fir]V pz ' 0.10) 


AS o dz 


and 


dB 




PZ _ 


AS o dz 


[« 


- I 8 it p Q to D a + D as dz . „ Q d2 


d5 o , , dD AS 


D AS 6 o 


^* 5 o ) 


)]b , 

dz o J pz 


(3.11) 


where 5 Q = d(B Q )/dz, and 


AS 


= D, 


2 2 

j 6 a . 


(3.12) 


Both equations (3.10) and (3.11) show two singular points where 
d AS = 0 or 6 q = 0; for incompressible medium ( 6 A -*• 0) the first 
singular point occurs when D A = 0 (see equation 3.1). This case 
is discussed by Tataronis and Grossmann (1973) who showed that 
the existence of singular points leads to a continuous spectrum 
of frequencies, in which case the normal mode analysis fails, and 
the solution is given by an integral over some range of wave 
frequencies in addition to a possible sum of discrete modes. 

They also showed that the continuous spectrum of frequencies is 
responsible for decaying of MHD waves by phase mixing (see also 
Lee and Roberts, 1986) . 
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V7e begin with the first singular point and calculate the 
phase velocity (U A = o»/k) for which the condition D AS = 0 is 
fulfilled. Introducing the angle $ between the magnetic field 
and the wave vector (see Figure 3)/ we obtain 


2 V S V A COS * 


ur = 


(3.13) 


V 


+ V 


Note that for <j> = 0 this expression describes the phase velocity 
of purely acoustic waves (or MHD slow waves when V A >> V s ) guided 
by an intense magnetic flux tube and is obtained when a thin flux 
tube approximation is considered (Defouw, 1976; Musielak, et al., 
1987) . 


We solve equations (3.10) and (3.11) at the point D AS =0 (a 
singular point solubility criterium) and find that for all layers 
in which |S q | goes uniformly through a minimum and shows no 
rotation there are no real wave frequencies for which the 
solutions can be obtained and only continuous spectrum solutions 
exist. In general, however, when rotation is taken into account 
the critical condition for existence of the solutions is found as 


„ 2 = -<2V 2 + V 2 U 2 - 2 (v| + V 2 U a X b tan*. 


(3.14) 


where X fl = d<j>/dz. This condition shows that non-zero and 
negative gradient of rotation of the magnetic field (X < 0) 

d 

introduces a possibility of finding a real wave frequency for 
which the condition (3.14) is satisfied (|x a tan<j>| > |x b |); then, 
we may define a real critical frequency o c (defined by the RHS of 
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equation 3.14) that satisfies D AS = 0 and normal mode solutions 
can be recovered . 

To consider the second singular point, we assume that the 
strength of magnetic field disappears at one particular point in 
the current sheet and that the first derivative of the square of 
the field strength with respect to z becomes zero in the vicinity 
of this point. The second derivative, however, remains non-zero 
leading to smooth and continous variations of B Q across the 
current sheet. In this case, 5 Q = 0 and using again singular 
point solubility criterium, we find 


2 

a) 


_1 

8 " p o 


dz * 


(3.15) 


If the second derivative of the magnetic field with respect to z 
is positive, then no real wave frequency can satisfy equation 
(3.15). The situation is different when the magnetic field 
exhibits a maximum inside the current sheet and the second 
derivative of the field becomes negative; it leads to the 
critical frequency o c defined by the RHS of equation (3.15) for 
which normal modes are the expected solutions of the wave 
equations . 


(c) Absence of Singular Points 

As described in the previous subsection the general MHD 
bending wave equations show two singular points which lead to 
limits on normal mode solutions. For an incompressible medium 
and non-rotating magnetic field, these two points are not 
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removable and exist for all wave frequencies that satisfy the 
condition D A = 0 in any layer (e.g. Tataronis and Grossmann, 

1973; Lee and Roberts, 1986). In the case considered in this 
paper, there are two agents; compressibility and rotation of the 
magnetic field vector which may remove singular points from the 
wave equations. 

Rotation of the magnetic field across the current sheet 
means that B Q (z) need not be zero somewhere in a "neutral" sheet 
and removes the second singular points from the wave equations. 
Observational data (Behannon et al., 1981) show that zero 
strength of the magnetic field apparently rarely exists for the 
HCS, mainly because of rotation, although variations of the 
magnetic field across the HCS often show a minimum. 

Now, we look for physical conditions necessary to remove the 
first singular point. Suppose that the condition (3.13) is not 
satisfied anywhere in the current sheet. Then, the phase 
velocity of transverse bending waves is equal to one of the local 
Alfven velocities characteristic of the current sheet and 
condition (3.13) is never fulfilled in the sheet, at least as 
long as the sound speed is less than the Alfven velocity (V A > 

V s ) . The latter inequality is not satisfied in the current sheet 
only in the vicinity of the point where the magnetic field 
disappears (i.e., a neutral sheet) or if the field variations 
across the sheet show a deep minimum. However, as shown above 
both latter cases can be prevented when the magnetic field 
rotates. It allows the condition V A > V g to be satisfied 
everywhere in the current sheet (see also Behannon et al., 1981) 
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and leads to absence of the first singular point. Note that when 
V A < V s the first singular point is not removable from the wave 
equations. Absence of the singular points leads to normal mode 
solutions for any wave frequency. 


(d) General Solutions 

Having discussed the singular points of the wave equations 
and their removal, we now present the general solutions of these 
equations valid for absence of the singular points (D AS * 0 
and S q * 0). To integrate the wave equations, we separate the 
variables and the physical parameters in equations (3.10) and 
(3.11), and obtain the solution for U pz in the form 




C C 
L AS o 

AS o 


exp 


f- « 2 f-^ 

L j V A *b°AS 



(3.16) 


the solution for B_„ is given by 

v* 


B = B 
pz pz 


C AS C o ^-Sq 1 

d as 6 o 


exp 


B 


[- » 2 [-s 5 *— <4 

L J V A VaS J 


(3.17) 


where U° , B° , C AO , C n and C are constant of integrations. The 

pZ pz Ao d O 

integrals can be evaluated explicitly only for very special cases 
that are discussed in the next subsection. 


(e) Simplifying the Problem: Limiting Cases 

In order to obtain some simple analytical solutions, we 
consider limiting cases for which both MHD bending waves 
equations can be significantly reduced. The simplest solutions 
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of equations (3.8) and (3.9) are found either when the gradient 
of magnetic field is neglected or when the dispersion relation 
for Alfven waves is satisfied (see two next cases). Simple 
solutions can also be obtained when the sound speed is much 
larger than the Alfven velocity (see third and fourth cases 
discussed below) . The problem would also be simplified if one 
could assume constant Alfven velocity across the current sheet; 
however, in the approach presented here (an isothermal medium) 
the pressure balance for the current sheet does not allow us to 
make this assumption. 

(i) rotates, B q (z) = constant 

In this simple case, both density and magnetic field 
gradients disappear ( X^ = X r = 0), leading to constant density 
across the current sheet. The rotation of the magnetic field 
does not introduce additional forces into the equation of motion 
as 


(V x 6 q ) x $ o = 0, 


(3.18) 


and the solutions for both Up Z and B pz are given by the 
dispersion relation for Alfven waves 


“ 2 ■ <*-V 2 = °* 


(3.19) 


which is valid for any perturbations of the velocity and magnetic 


field 
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Purely transverse Alfven waves are the only MHD bending 
waves for a current sheet consisting of an arbitrary rotation of 
the magnetic field through any angle, including 180° for a 
neutral sheet. To calculate the phase velocity U A for these 
waves, we introduce the angle 0 between the wave vector It and the 
y-axis (the axis is fixed in the direction of the magnetic field 
at the upper edge of the current sheet) and obtain 


U A = ± V A cos ( o( z) - 9) . 


(3.20) 


For propagation along the y-axis (0 = 0°) and arbitrary direction 
of the magnetic field with respect to the y-axis one gets 
U A = V A cos a; and, for a magnetic field parallel to the y-axis 
and arbitrary direction of the wave propagation equation (3.21) 
gives U A = V A cos 0. In addition, one may introduce the 
angle <j> between the magnetic field and wave-vector ( 4> = a - 0 ) to 
further simplify the expression for the phase velocity. All 
latter results are well-known solutions for Alfven waves in a 
homogeneous medium with uniform magnetic field. 


(ii) S Q rotates, D A = 0 

In this very special case, the dispersion relation for 
transverse Alfven waves is satisfied and the solution for U 0 _ 
(equation 3.16) reduces to the form 


U 


pz 


= 


C C V 
AS o S 


pz 


2 B' 


2 2 
V V A 


(3.21) 
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This solution is valid only for one particular point in the 
current sheet where D A = 0 for each wave frequency; note that 
this point is not a singular point for the wave equations (3.8 
and 3.9) and normal mode solutions exist. In general the Alfven 
velocity varies across the sheet and the dispersion relation is 
not satisfied. This means that the global solution (3.16) for 
U pz must equal this value at the critical layer where D A = 0 . 

(iii) rotates, B Q (z) is arbitrary, V s >> V A 

The later assumption { 0 A 0 ) is equivalent to neglecting 
the density gradient in the continuity equation (see Bertin and 
Coppi 1985) and leads to the incompressible fluid approximation 
(V«t) = 0). Formally, this Boussinesq type of assumption makes 
the approach of studying MHD bending waves inconsistent, mainly 
because the considered medium is compressible and the density 
gradient in the continuity equation cannot be neglected. 
Referring to the results previously obtained, we discuss this 
case in order to present the solution of the wave equation for 
U pz given by equation (3.16). For S A -*• 0 the behaviour of MHD 
bending waves is described either by the wave equation (3.1) or 
by equation (3.16) reduced by the assumption D AS = D A . Now, the 
first critical point is defined by the condition D A = 0 and 
cannot be removed by the rotation alone . This critical point is 
always present in equation (3.16), reduced for an incompressible 
medium, and an integral over some range of wave frequencies 
appears in addition to a sum of discrete modes (Tataronis and 
Grossman, 1973; Priest, 1982). Simple solutions are found only 
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for real wave frequencies that do not satisfy D A = 0 anywhere in 
the current sheet (however the condition can be satisfied outside 
the sheet) and according to equation (3.16), we obtain 


U 


pz 



c c 

A ^O 


2 Vb°A 



(3.22) 


where C A is a constant of integration? this integral can be 
calculated when simple variations of the magnetic field across 
the current sheet are assumed (see the next case). 


(iv) rotates, B Q (z) = B e tanh (z/d) and Vg » V A 

In order to evaluate the integral in equation (3.24), we 
have to calculate the Alfven velocity as well as x^ for this 
specific variation of the magnetic field (Figure 4). From the 
pressure balance, and assuming 2 ^ = B and 

G G G 

= pf 2 ^ = P we obtain 
e e e 
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(3.23) 


where 8 Ae = v Ae / v g» v Ae = B e /^np e ) ^ and 2d is a thickness of 
the current sheet. It can also be shown that 


(P Q " Pe } cosh 2 (z/d) = I p e s 2 e 


= const. 


(3.24) 


Using equation (3.25), the Alfven velocity is found in the form 


v: = 2 v 


sinh^ ( z/d ) 


Ae 8 2 + 
s Ae 


2 cosh (z/d) 


(3.25) 
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which can be further simplified by the assumption V g >> 

V A (8 Ae -*■ O'). Then, a simple analytical formula is found for the 
integral in equation (3.22) and the solution is given by 



yO C A^o cosh ^ ( z/d ) 
pz 2D A d s £ n h K ( 2 /(j) 


coth( z/d) 


3 

exp[- V 2 < cosh -3 (z/d)] 

(3.26) 


2 2 2 

where d Q is a constant of integration and k = <d /V A d . The 
solution given by equation (3.26) is valid only for those 
frequencies which do not satisfy the dispersion relation D A = 0 
anywhere in the current sheet; otherwise, the continuous spectrum 
of frequencies is the expected solution of the problem. 


4. DISCUSSION 

We have presented the normal mode analysis for ideal MHD, 
incompressible (the Alfven type of purely transverse) bending 
waves on a two-dimensional current sheet. The current sheet is 
of completely general structure, in equilibrium with its 
homogeneous surroundings. The magnetic field may undergo both 
amplitude and directional changes within the sheet to simulate 
what is observed to occur in the heliospheric current sheet 
( HCS ) . The behavior of MHD incompressible waves is fully 
described by the wave equations (3.8) and (3.9) which are the 
most important results of this paper. In the limit of no field 


rotation and infinite sound speed, we recover the results of Lee 
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and Roberts (1986) and others; that the normal mode analysis 
fails, there is a continuous spectrum of solutions for any given 
wave number, and that waves undergo decay through "phase mixing" 
or "mode conversion" - ultimately involving dissipative and 
nonlinear processes that are beyond the scope of our work but 
which have analyzed by others ( Ste inolfson , et al., 1986, op 
ci t . ) . 

A new result is that in the presence of field rotation, 
finite sound speed, and/or a local maximum. in field strength, it 
is possible to recover conditions under which normal modes 
exist. Furthermore, there are specific conditions under which 
the singularity in the wave equation that leads to continuous 
spectrum solutions can be removed by using a critical point 
solubility criterium such that a different class of normal modes 
is recovered. There are even conditions under which continuous 
spectrum solutions disappear altogether. 

These results mean that disturbances on a current sheet need 
not always decay through mode conversion - which leads us to 
speculate on the physical processes in the HCS. First, we note 
that the HCS at 1 AU apparently generally attains a field 
reversal through field rotation and hence may not support decay 
through mode conversion. Second, "ripples" are commonly observed 
on the HCS. These ripples are, therefore, probably Alfven-like 
normal modes which are undergoing no further decay. However, 
near the sun the establishment of the HCS does not necessarily 
require a field rotation. Thus waves present near the sun - and 
it is reasonable to assume waves are present in that they are 
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generally observed at 1 AU - would undergo mode conversion and 
decay, ultimately leading to dissipation and local magnetic 
diffusion such that the magnetic field relaxes into a 
configuration having a field rotation and the waves would no 
longer decay and diffusion ceases. However, this scenario is not 
completely in agreement with the data for a combination of 
reasons. First, Steinolfson, et a 1. (1986) report that viscous 
decay is more efficient than ohmic decay under coronal 
conditions. In view of this, we would expect heating in the 
vicinity of the HCS. However, Borinni, et al . (1981) report just 
the opposite - the proton temperature tends to be a minimum in 
the vicinity of the HCS. The only way to circumvent this 
discrepancy is if dissipation in the interplanetary medium 
behaves differently than in the corona so that most of it occurs 
via a decay process involving reconnection of the magnetic field. 

Returning to the analysis of Bertin and Coppi (1985), it is 
clear that their idea of MHD bending waves on the HCS now 
deserves considerable attention, even though they misapplied 
their idea in attempting to explain the warping of the HCS. Data 
show that warping of the HCS can be explained as the imprint of 
coronal structure. However, bending waves analyzed for small- 
scale waves, when the WKB approximation universally used in their 
anlaysis is valid, offer the possibility of understanding the 
character and evolution of the ripples known to exist on the 
HCS. Furthermore, our analysis, which includes all types of 
configurations describable by a two-dimensional layer whether or 
not the field actually reverses, offers tools for defining when 
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surface waves will decay via mode conversion in the general 
astrophysical context. Finally, the new closed integral solution 
we have found is an analytic tool that simplifies finding 
solutions under those conditions for which normal modes exist. 
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FIGURE CAPTURES 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Topology of the heliospheric current sheet using 
Wilcox Solar Observatory line-of-sight magnetic field 
data, a potential field model of the corona with a 
source surface radius of 2.5 R Q , and a 400 km/s solar 
wind speed. At position A, Halley's Comet intersects 
the current sheet on 13.1 March 1986 - at the same 
time and place as the ESA Giotto spacecraft. The 
Earth was at position B at the same time. 

Sketch of the model layer with definition of 
coordinate system (2a) and variations of the magnetic 
field in the current sheet plane (2b). 

. Definition of the pertinent symbols relating the 
directions of the local magnetic field and the wave 
vector with respect to the coordinate system. 

Variations of the magnetic field $ (z) = § tanh (z/d) 

U C 

in the current sheet plane; 2d is a thickness of the 
current sheet. The semi-circular dashed lines show 
the path of the tip of the & Q (z) vector with two 
arbitrary intermediate vectors shown 
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